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r-». ! We calculate the nuclear cross section for vector meson exclusive photoproduction within 

the QCD color dipole picture and in the Regge approach. For the former approach, we have 
considered the phenomenological saturation model, whereas for the latter we use a model based 
on the dipole Pomeron framework. Theoretical estimates for scattering on both light and heavy 
nuclei are given over a large range on energy. 
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1 Introduction 



Exclusive vector meson production by real and virtual photons is an outstanding process pro- 
viding important information on the transition region from the soft dynamics (at low virtualities 
of the photon Q 2 ) to the hard perturbative regime at high Q 2 [1, 2]. In principle, a perturbative 
approach is only justified if a hard scale is present in the process, e.g. the photon virtuality 
and/or a large mass of the vector meson. For photoproduction of light mesons, such scale is 
not present and one has to rely on non-perturbative models. In general, a simple Regge pole 
phenomenology, with a soft Pomeron having intercept larger than one, is enough to describe 
the energy dependence of the meson cross section at the present accelerators. On some pQCD 
approaches, as the saturation model, even this soft process can be described, where the transi- 
tion is set by the saturation scale. Both models give an effective Pomeron intercept increasing 
with photon virtuality and meson mass. For <fi and J/^> mesons, the Pomeron intercept is 
considerably large and consistent with that one obtained in usual pQCD approaches. Despite 
the good agreement for the currently available energies, an extrapolation to higher energies of 
the experimental fits implies a large growth for the cross section which it would violate the 
unitarity at sufficiently high energies. Therefore, dynamical modifications associated to the 
unitarity corrections are also expected to be present in this process [3, 4]. Moreover, these 
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effects should be enhanced in nuclear processes [5, 6]. In particular, in the planned eA colliders 
at HERA and RHIC, the experimental analyzes of the exclusive vector meson production could 
be very useful to constrain the QCD dynamics [7, 8] . 

Our goal in this paper is to investigate the high energy vector meson exclusive photopro- 
duction on nuclei. In particular, we improve the previous analyzes of vector meson production 
[9] which are based on an extrapolation of the DESY-HERA experimental fits for the proton 
case. Here, we will consider unitarized cross sections ab initio. In order to do so, we con- 
sider two distinct and well established theoretical scenarios, which do not violate the unitarity 
bound in the asymptotic regime to be probed in future colliders. This allows us to analyze the 
nuclear vector meson protoproduction as a potential process to discriminate between these dif- 
ferent theoretical approaches. First, we consider the color dipole description of the 7 A — > V A 
(V = p, ou, 0, J/^) process, which is quite successful for the proton case [3, 4] and can be 
extended to nuclei targets via Glauber-Gribov formalism. It is important to quote the pio- 
neering papers [10, 11, 12, 13, 14], where these issues were first addressed and their further 
developments [15, 16, 17] as well. Related calculations in the /cj_-factorization approach (dipole 
approach is equivalent to it at leaging logarithmic aproximation) can be found in the Refs. 
[18, 19, 20, 21]. In the color dipole approach, the degrees of freedom are the photon (color 
dipole) and meson wavefunctions as well as the dipole-nuclei cross section. Such an approach 
enables to include nuclear effects and parton saturation phenomenon. The latter one is char- 
acterized by a typical momentum scale Q SSut (saturation scale) and it has been constrained by 
experimental results in deep inelastic scattering (DIS) and diffractive DIS [22]. Here, we will 
use an extension of the phenomenological saturation model for nuclear targets [23] . This model 
reasonably describes the experimental data for the nuclear structure function and has been used 
to predict the nuclear inclusive and diffractive cross sections for heavy quark photoproduction 
[24]. The nuclear saturation scale, Q s a, provides the transition between the color transparency 
and the saturation regimes in the nuclear scattering. Concerning vector meson production, 
our starting point is the recent work in Ref. [25], where different meson wavefunctions and 
dipole cross sections are considered for the proton case. It is worth mentioning that although 
light meson photoproduction to be a soft process by definition, it is consistently described in 
the QCD color dipole picture whether there is a suitable model for the soft-hard transition, as 
occurring in the the saturation model. 

In order to compare the saturation approach with a successful nonperturbative formalism, 
we consider a Regge inspired model given by the dipole Pomeron framework [26]. The reason 
for this particular choice is that in this model the soft Pomeron having intercept equal one, thus 
it does not violate unitarity for hadron-hadron and vector meson production at higher energies. 
Moreover, it describes with good agreement hadronic cross section and even DIS data in a wide 
range of photon virtualities [27]. For meson production on proton target, our starting point 
is the recent work of Ref. [28], where exclusive photoproduction by real and virtual photons 
is described with good agreement. The extension to nuclei is provided by the assumption of 
vector meson dominance and the Glauber-Gribov formalism. 

This paper is organized as follows. In the next section, we present a brief review of the 
exclusive meson production in the color dipole picture for proton target and its extension 
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to the nuclear case. For the dipole-proton (nucleus) cross section we have considered the 
phenomenological saturation model, which is shortly described. In Sec. 3, the dipole Pomeron 
parameterization for vector meson photoproduction is presented and its extension to nuclear 
targets is considered. The results coming out from both models are presented and discussed in 
Sec. 4. Finally, in Section 5 we summarize our conclusions. 



2 Vector meson production in the color dipole approach 

Let us introduce the main formulas concerning the vector meson production in the color dipole 
picture. First, we consider the scattering process •yp — > Vp, where V stands for both light and 
heavy mesons. Further, one extends this approach to the nuclear case. The scattering process 
can be seen in the target rest frame as a succession in time of three factorizable subprocesses: 
i) the photon fluctuates in a quark-antiquark pair (the dipole), ii) this color dipole interacts 
with the target and, iii) the pair converts into vector meson final state. Using as kinematic 
variables the 7*iV c.m.s. energy squared s = W 2 N — (p + q) 2 , where p and q are the target and 
the photon momenta, respectively, the photon virtuality squared Q 2 = —q 2 and the Bjorken 
variable x = Q 2 /(W 2 N + Q 2 ), the corresponding imaginary part of the amplitude at zero 
momentum transfer reads as [29], 

Xm A (tp ^Vp) = Y._jdz d 2 r ^(z, r, Q 2 ) a%^{x, r) ^(z, r) , (1) 

h,h 

where ^ hh (z, r) and ^/^(z, r) are the light-cone wavefunctions of the photon and vector me- 
son, respectively. The quark and antiquark helicities are labeled by h and h and reference to 
the meson and photon helicities are implicitly understood. The variable r defines the relative 
transverse separation of the pair (dipole) and z (1 — z) is the longitudinal momentum frac- 
tions of the quark (antiquark). The basic blocks are the photon wavefunction, \1/ 7 , the meson 
wavefunction, ^t,l-> an d the dipole-target cross section, o"^ get . 

In the dipole formalism, the light-cone wavefunctions ^ h ^(z, r) in the mixed representation 
(r, z) are obtained through two dimensional Fourier transform of the momentum space light- 
cone wavefunctions '^ h - h (z ) k) (see more details, e.g. in Refs. [1, 4, 25]). The normalized light- 
cone wavefunctions for longitudinally (L) and transversely (T) polarized photons are given by 
[30]: 

Kn^r) = f-^5^_- h ee f 2z(l-z)Q^l, (2) 



*^\*>r) = ±^ee / [ie ± ^(A ±) ^-(l-z)^J9 r + m / VM 1 ^^ ( 3 ) 

where e 2 = z(l — z)Q 2 + m 2 . The quark mass rrif plays a role of a regulator when the 
photoproduction regime is reached. Namely, it prevents non-zero argument for the modified 
Bessel functions K 0t i(er) towards Q 2 — * 0. The electric charge of the quark of flavor / is given 
by ee f . 




K (er) 



3 



For vector mesons, the light-cone wavefunctions are not known in a systematic way and 
they are thus obtained through models. The simplest approach assumes a same vector current 
as in the photon case, but introducing an additional vertex factor. Moreover, in general a 
same functional form for the scalar part of the meson light-cone wavefunction is chosen. Here, 
we follows the analytically simple DGKP approach [30], which is found to describe in good 
agreement vector meson production as pointed out in Ref. [25]. In this particular approach, one 
assumes that the dependencies on r and z of the wavefunction are factorised, with a Gaussian 
dependence on r. The DGKP longitudinal and transverse meson light-cone wavefunctions are 
given by [30], 



z(l- z) 5 h - h f L (z) exp 



u 2 L r 2 



± 



2y/N c e f 

[zS h± - hT - (1 - z)5 hTrh± ] + — S h± - h± 



(4) 



x^==—f T (z) exp 



-uj 2 L r 2 



(5) 



V2N~ c e f 

where if is the effective charge arising from the sum over quark flavors in the meson of mass my. 
The following values if = l/y/2, l/3\/2, 1/3 and 2/3 stand for the p, u, <f) and J/^, respectively. 
The coupling of the meson to electromagnetic current is labeled by fy = 3 m v Y e + e - /A na 
(see Table 1). The function /t,l(z) is given by the Bauer-Stech-Wirbel model [31]: 

m 2 v (z-1/2) 2 



2 

em 



fr,L(z) = Nt,l \jz{l - z) exp 



2oo 2 



T, L 



(6) 



The meson wavefunctions are constrained by the normalization condition, which contains 
the hypothesis of a meson composed only of quark- ant iquark pairs, and by the electronic decay 
width ry^ e + e -. Both conditions are respectively given by [32, 4], 



Ejd 2 rdz\^\z, r)\ 2 =l, 



(7) 



h,h 



5/ 



d 2 r 



dz 



h,h 



(2tt) 2 z(1 -z) 



(1 - z)Q 2 + k 2 + m)\ *l- h (k, z)*l*- h (k, z) = efymv « • e v ) . (8) 



The above constraints when used on the DGKP wavefunction produce the following relations 

[25], 



f 

JO 



^2W c e } 

dz z(l - z) f L (z) = [ 
Jo 



i 2[z 2 + (1- z) 2 ]ul + m) 



dz 



2m\r z{l — z) 



h{z) 



where 



h = 



dzz 2 (l-zYfi(z) 



-I 



1 , ^ + (1 

dz 



uj\ + m 2 2 



m 2 v 



m ■ 



(9) 
(10) 

(11) 
(12) 
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V(m v ) 


ev 


/v 


UJt 


Mr 


By 


MeV 




[GeV] 


[GeV] 




[GeV- 2 ] 


p(770) 


1/V2 


0.153 


0.218 


8.682 


9.00 


u (782) 


1/3^2 


0.0458 


0.210 


10.050 


10.14 


0(1019) 


1/3 


0.079 


0.262 


8.000 


8.92 


J/V (3097) 


2/3 


0.270 


0.546 


7.665 


4.57 



Table 1: Parameters and normalization of the DGKP vector meson light-cone wavef unctions. 
Results obtained using quark mass values from the saturation model (see text). 

The relations in Eq. (9) come from the normalization condition, whereas the relations in 
Eq. (10) are a consequence of the leptonic decay width constraints. The parameters ut,l and 
jVt,l are determined by solving (9) and (10) simultaneously. In Tab. 1 we quote the results for 
the transverse component, which it will be used in our further analysis in the photoproduction 
case (longitudinal component does not contribute at Q 2 = 0). To be consistent with the 
phenomenological saturation model, which we will discuss further, we have used the quark 
masses m u ^ s = 0.14 GeV and m c = 1.5 GeV. In the case of meson, we follow Ref. [30] and 
take m s = m u ^ d + 0.15 GeV. We quote Ref. [25] for more details in the present approach and 
its comparison with data for both photo and electroproduction of light mesons. 

Finally, the imaginary part of the forward amplitude can be obtained by putting the expres- 
sions for photon and vector meson (DGKP) wavef unctions, Eqs. (2-3) and (4-5), into Eq. (1). 
Moreover, summation over the quark/antiquark helicities and an average over the transverse 
polarization states of the photon should be taken into account. The longitudinal and transverse 
components are then written as [25, 30] 

XmA L = J d 2 r dz y /a~f v 2z 2 ( y l-z) 2 f L (z) exp 

d 2 r dzy/a^ f v f T (z) exp — | — 

x {^V + (1 - z) 2 } Kl (er) + ^K^er)) r) , (14) 

[ my my J 

with <r^p get being the dipole-proton cross section in the nucleon case and the dipole-nucleus 
cross section for scattering on nuclei. For the proton case, there are a lot of phenomenology 
for p and J/ty production using recent pQCD parameterizations for the dipole-proton cross 
section [25] or considering nonperturbative QCD calculations based on stochastic QCD vacuum 
[33, 34]. In the next subsection, we briefly review the dipole-nucleon (nucleus) given by the 
phenomenological saturation model, which it will be considered in our numerical estimates. 

In order to obtain the total cross section, we assume an exponential parameterization for 
the small \t\ behavior of the amplitude. After integration over |t|, the total cross section for 



Q^oM^ ct (x,r),(13) 
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vector meson production by real/virtual photons in the nucleon (proton) case reads as, 



a (jp — > Vp) 



[ImA(s, t = 0)} 2 
16irB v 



(1 + /3 2 ) 



(15) 



where (3 is the ratio of real to imaginary part of the amplitude and By labels the slope parameter 
(we quote the values we have used in Table 1). The values considered for the slope parameter 
are taken from the parameterization used in Ref . [3] . For the p case, we have taken a different 
value in order to describe simultaneously HI and ZEUS photoproduction data. 

In addition, Eqs. (13-14) represent only the leading imaginary part of the positive-signature 
amplitude, and its real part can be restored using dispersion relations Tie A = tan(7rA/2) XmA. 
Thus, for the (5 parameter we have used the simple ansatz, 



with A c ff = A e ff (W 7 jv, Q 2 ) the effective power of the imaginary amplitude, which depends on 
both energy and photon virtuality. The correction coming from real part in photoproduction, 
where only transverse component contributes, is about 3% for light mesons and it reaches 
13% for J/\l/ at high energies. It is worth mentioning that a different computation of the (3 
parameter, as in Ref. [25], produces a larger effect even in the photoproduction case. An 
additional correction is still required for heavy mesons, like J/^> . Namely, skewedness effects 
which takes into account the off-forward features of the process (different transverse momenta 
of the exchanged gluons in the t-channel), are increasingly important in this case. Here, we 
follow the studies in Ref. [35], where the ratio of off-forward to forward gluon distributions 
reads as [35], 



In the case of nuclear targets, By is dominated by the nuclear size, with B ~ R\ (Ra = 
1.2 A 1 / 3 fm is the nuclear radius) and the non- forward differential cross section is dominated 
by the nuclear form factor, which is the Fourier transform of the nuclear density profile. Here 
we use the analytical approximation of the Woods-Saxon distribution as a hard sphere, with 
radius Ra, convoluted with a Yukawa potential with range a = 0.7 fm. Thus, the nuclear form 
factor reads as [9] , 





d ln[ImA(s, t = 0)] 
d Ins 



(16) 




(17) 



F(q 




47rpo 
Aq 3 



[sm(qR A ) - qR a cos(qR A )} 



1 + a 2 q 2 



1 



(18) 



where p = 0.16 fm 3 . 

The photonuclear cross section is given by 



a( 1 A^VA)= [ImA ^ t = 0)]2 



(1 + /5 2 ) / dt\F(t)\ 2 , 



(19) 
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with t min = (my/2u) 2 , where u is the photon energy. 

Having introduced the main expressions for computing vector meson production in the 
color dipole approach, in what follows we present the saturation model and its extension for 
the scattering on nuclei targets. 



2.1 Dipole- nucleus cross section in the saturation model 

For electron-proton interactions, the dipole cross section cr^ r ° ton , describing the dipole-proton 
interaction, is substantially affected by nonperturbative content. There are several phenomeno- 
logical implementations for this quantity. The main feature of these approaches is to be able to 
match the soft (low Q 2 ) and hard (large Q 2 ) regimes in an unified way. In the present work, we 
follow the quite successful saturation model [22], which interpolates between the small and large 
dipole configurations, providing color transparency behavior, a^ v ~ r 2 , as r — > and constant 
behavior, a dip ~ cr , at large dipole separations. The parameters of the model have been ob- 
tained from an adjustment to small x HERA data. Its parameter- free application to diffractive 
DIS has been also quite successful [22] as well as its extension to virtual Compton scattering 
[36], vector meson production [3, 25] and two-photon collisions [37]. The parameterization for 
the dipole cross section takes the eikonal-like form [22], 



proton/- 2\ _ 
a dip \ X i r ) - °0 



= GeV 2 , (20) 



where the saturation scale Q 2 at defines the onset of the saturation phenomenon, which de- 
pends on energy. The parameters, obtained from a fit to the small- x HERA data, are a = 
23.03 (29.12) mb, A = 0.288 (0.277) and x = 3.04 • 10~ 4 (0.41 • 10" 4 ) for a 3-flavor (4-flavor) 
analysis. An additional parameter is the effective light quark mass, rrif = 0.14 GeV, which 
plays the role of a regulator for the photoproduction (Q 2 = 0) cross section, as discussed be- 
fore. The charm quark mass is considered to be m c = 1.5 GeV. A smooth transition to the 
photoproduction limit is obtained via the scaling variable [22], 

Q 2 + 4:m 2 

i= Q^w£- (21) 

The saturation model is suitable in the region below x = 0.01 and the large x limit needs 
still a consistent treatment. Making use of the dimensional-cutting rules, here we supplement 
the dipole cross section, Eq. (20), with a threshold factor (1 — x) nthICS , taking n tn res — 5 for a 
3-flavor analysis and n tri rcs = 7 for a 4-flavor one. This procedure ensures consistent description 
of heavy quark production at the fixed target data [38] . 

Let us discuss the extension of the saturation model for the photon-nucleus interactions. 
Here, we follow the simple procedure proposed in Ref. [23], which consists of an extension to 
nuclei of the saturation model discussed above, using the Glauber-Gribov picture [39], without 
any new parameter. In this approach, the nuclear version is obtained replacing the dipole- 
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nucleoli cross section in Eq. (1) by the nuclear one, 



.nucleus 

dip 



(x, r 2 ; A) = 2 J ' d 2 b j 




.proton 
dip 



(x, r 2 )] } 



(22) 



where b is the impact parameter of the center of the dipole relative to the center of the nucleus 
and the integrand gives the total dipole-nucleus cross section for a fixed impact parameter. 
The nuclear profile function is labeled by Tjy(b), which will be obtained from a 3-parameter 
Fermi distribution for the nuclear density [40]. The above equation sums up all the multiple 
elastic rescattering diagrams of the qq pair and is justified for large coherence length, where 
the transverse separation r of partons in the multiparton Fock state of the photon becomes 
as good a conserved quantity as the angular momentum, i. e. the size of the pair r becomes 
eigenvalue of the scattering matrix. It is important to emphasize that for very small values of 
x, other diagrams beyond the multiple Pomeron exchange considered here should contribute 
(e.g. Pomeron loops) and a more general approach for the high density (saturation) regime 
must be considered. However, we believe that this approach allows us to obtain lower limits 
of the high density effects at eRHIC and HERA-A. Therefore, at first glance, the region of 
applicability of this model should be at small values of x, i.e. large coherence length, and 
for not too high values of virtualities, where the implementation of the DGLAP evolution 
should be required. Consequently, the approach is quite suitable for the analysis of exclusive 
vector meson photoproduction in the kinematical range of the planned lepton-nucleus colliders 
(eRHIC and HERA-A). Furthermore, it should be noticed that the energy dependence of the 
cross sections is strongly connected with the semi-hard scale (the saturation momentum scale). 
Namely, the saturation effects are larger whether the momentum scale is of order or larger than 
the correspondent size of the vector meson and the energy growth of the cross section is then 
slowed down. 

3 Vector meson photoproduction in the dipole Pomeron 
framework 

Let us summarize the main features and expressions for the nonperturbative approach given 
by the dipole Pomeron model [26]. This model describes the vector meson exclusive photopro- 
duction [28] data from HERA without need of a Pomeron contribution with intercept higher 
than 1, thus not violating the Froissart-Martin bound. The picture of the interaction is given 
by a photon fluctuating into a quark- ant iquark pair and further the nucleon (proton) interacts 
with it through Pomeron or secondary Reggeon exchange. After that, the pair converts into a 
vector meson. In general lines, this picture is quite similar to that one for interaction among 
hadrons in the Regge limit. In particular for photoproduction, the representation of the photon 
as a hadron is reasonably supported and the Regge pole theory, with a Pomeron universal in 
all hadron-hadron interactions, can be safely used there. 

For the Pomeron contribution we follow Ref. [28] and one considers the dipole Pomeron, 
which gives a very good description of all hadron-hadron total cross sections. As Pomeron 
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and secondary Reggeons would be universal objects in Regge theory, the corresponding j- 
singularities of photon-proton amplitudes and their trajectories at the photoproduction limit 
coincide with those appearing in pure hadronic amplitudes. In particular, the restriction on the 
Pomeron intercept implied by the Froissart-Martin bound suggests that it is a more complicated 
singularity instead of a simple pole having a universal intercept otp > 1, which one would apply 
also to DIS. In the case of a dipole Pomeron, it is a double j-pole leading to cx^' 7/l oc Ins and 
unitarity requirements are covered. 

Let's consider the usual Mandelstam variables, s = W 2 N = m 2 N + 2m N u — Q 2 {m N is 

the nucleon mass), t (momentum transfer) and defining the scaling variable, Q = Q 2 + m v . 
Making use of the latter quantity, the dipole Pomeron model can be generalized for virtual 
external particles. The scattering amplitude is given by the contribution of Reggeons at low 
energies and the dipole Pomeron (ajp(t = 0) = 1) dominates at higher energies. A simple pole 
parameterization is used for the /-Reggeon. The parameters ap(t) and ant(t) are universal 
and do not depend of the reaction, whereas couplings gi, energy scales soi and slopes hi are 
functions of the scaling variable and the same for all reactions. 

Taking into account the features discussed above, the differential elastic cross section is 
written as [28] 

(23) 

where in the photoproduction case the amplitudes for the secondary Reggeons and Pomeron 
contributions are parameterized as [28], 

Ajp(W 2 ,t; m 2 v ) 



A M (W 2 N ,t] m 2 v ) 

where one takes a linear Pomeron trajectory ajp(t) = 1 + «jp(0) t, with the usual value for the 
slope a'jp(0) = 0.25 GeV~ 2 . The Reggeons and Pomeron couplings are written as, 

#o,i (t; m v ) 

9m (t ; m v ) 

where the couplings (? ,i, the energy scale Wq (GeV 2 ) and t-slope bjp (GeV -2 ) for the Pomeron 
are adjustable parameters of the model; in addition m 2 is the proton mass. The notation V 
stands for p, <fi and J/ip, whereas M = /, 7r for u>. The remaining constants for the Reggeons, 
Qfi 9n, b 2 jR (GeV -2 ) are also adjustable parameters. One uses the same slope 6|j for / and n 



rli 



( 7 p _^ Vp) = 4vr A F (s,t 



m v ) 



A M (s,t] m 2 y) 



ig (t;my) y-i 
+ igi(t\m v ) In 
ig M (t] my) - 



W* N 



mz 



a P (t)-l 



W 2 



m 



v 



W* N -m 2 p 



W* N -m 2 p 



Wo + m v J V w o+ ™ 2 v 



m. 



w 2 lN 

W2 + ml 



2 \ a M (t)-l 



V 



a P (t)-l 



(24) 
(25) 



9o,imp 
(W 2 + m 2 v ) 2 

gnm 2 p 
(W 2 + m 2 ) 2 



cxp 



(bit 



exp (b^t 



(26) 
(27) 
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Reggeon exchanges. We quote Ref. [28] for details on the fit procedure and the tables for the 
fitted parameters. Let us stress that the only variable that differentiates among the various 
vector meson elastic cross sections is the mass of the vector mesons. 

Following Ref. [28], some additional comments are in order. The behavior in the thresh- 
old region is given by multiplying the amplitudes by the correction factor (1 — x) s , where 
x = (m p + m v ) 2 /W 2 N and (m p + my) is the reaction threshold. The power 5 = \jm Y /m\ 
drives the energy dependence in that region, with m\ (GeV 2 ) fitted from data. Concerning 
details in Ref. [28] when taking into account the secondary Reggeons, for p, and J /if) me- 
son photoproduction the scattering amplitude was written as the sum of a Pomeron and / 
contribution. Although according to Okubo-Zweig rule, the / meson contribution ought to be 
suppressed in the production of and J /if) mesons, the / meson contribution was added even 
in the J /if) meson case. For J /if), it was found to be negligible whereas it is sizable for meson 
production. 

The 7p — > Vp process can be used as input in the calculations of the total cross section for 
the reaction 7 A — > VA. A major simplification comes from the use of vector meson dominance, 
which allows to relate this photoproduction cross section to the cross section for the forward 
elastic Vp — > Vp scattering. Following vector meson dominance [41, 42], 



da(^p — > Vp) 



dt 



t=o 



47ra em d a(Vp — > Vp) 

~Tv " 



dt 



t=o 



(28) 



where, t is the squared 4-momentum transfer between the proton and vector meson, a em is 
the electromagnetic coupling constant and fy is the vector meson-photon coupling, fy = 
4n myal m / (3T v ^ e + e -), with my the vector meson mass and T v ^ e + e - the leptonic decay par- 
tial width. Values for fy/^it are given in Table II of Ref. [9]. Using the Optical theorem, the 
total cross section is given by, 



cr to t(Vp -> Vp) 



16tt 



da(Vp — > Vp) 



dt 



t=o 



(29) 



The scattering cross section from heavy nuclei can be found by a (quantum mechanical) 
Glauber-Gribov calculation, 



a tot (VA ^VA) = 2j d 2 b [l - exp (-^ T A (b) a tot (Vp 



Vp) 



(30) 



As referred before, the nuclear profile function is labeled by T^(b), which will be obtained from 
a 3-parameter Fermi distribution for the nuclear density [40]. The Optical theorem for nucleus 
A and vector meson dominance are then used to find the following relation, 



da^A -> VA) 



dt 



a 



t=o 



4/ 



em 2 

/T 



a 2 tot {VA - VA) . 



(31) 



From this equation one can directly understand the A-dependence in two limiting cases: in the 
transparent limit there is a A 2 behavior (typical of coherent processes) and in the black disc 
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Figure 1: The total cross section for p and u photoproduction on proton as well as for light (Ca) 
and heavy (Pb) nuclei. The solid lines stand for the QCD color dipole approach and the dashed 
ones for the soft dipole Pomeron approach. Experimental high energy data from DESY-HERA 
collider on proton target are also shown. 



limit we have an As rise with the nuclear number A. The total photonuclear cross section is 
then given by 



o^A - VA) 



daft A -> VA) 



dt 



t=o Jt 



dt\F(t)f 



(32) 



where F(t) is given in Eq. (18). As F(t) is A-dependent we have that integration over t yields 
a factor of A~s, which implies a As (As) behavior in the transparent (black disc) limit. 



4 Nuclear vector meson exclusive photoproduction 

In this section we compute the nuclear cross section for the exclusive photoproduction of vector 
mesons. Here, we compare the QCD approach given by the saturation model extended to 
nuclei targets within the color dipole picture as well as a nonperturbative approach rendered 
by the dipole Pomeron model, which does also not violate unitarity at high energies. We focus 
on the energy range and nuclei targets expected for the future lepton-nuclei colliders (eRHIC 
and eHERA) and also for available range to be covered in ultraperipheral heavy-ion collisions 
(UPC's) at LHC. 

In Figs. 1 and 2 are shown the results for the p, u), <fi and J/^ photoproduction cross 
section as a function of energy for different nuclei, including the proton case. The results in 
the pQCD and Regge approaches present mild growth on W^n at high energies stemming from 
the high energy behavior of the models, whereas the low energy region is consistently described 
through the threshold factor. For the proton, the experimental data from DESY-HERA collider 
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[43, 44, 45, 46, 47, 48] are also included for sake of comparison. We can see that the p, u and 
results vary only slowly with energy, in contrast with the J/^f predictions. The saturation 
model (solid lines) gives a flatter energy dependence in comparison with the dipole Pomeron 
model (dashed lines). On the other hand, for J/ty the situation changes, where the dipole 
Pomeron model producing a mild increasing at high energies in comparison with the saturation 
model. 

Let us discuss the results coming from the color dipole approach. For the proton case we 
have used Eqs. (14-15) and the dipole-proton cross section given by Eq. (20). The contribution 
of the real part of amplitude is small for the light mesons, whereas is sizeable for the J/\l/ case. 
Moreover, the skewedness correction to the J/\l/ photoproduction is important, providing a 
larger overall normalization, as discussed in Sec. 2. These features remain in the nucleus case, 
where we have used Eq. (19) and the dipole-nucleus cross section given by Eq. (22). It is 
worth mentioning that the effective power of the imaginary part of amplitude is slowed down 
in the nuclear case and this has implications in the corrections of real part and skewedness. 
The results for photonuclear production on nuclei is consistent with the studies in Ref. [9], 
except for J/^> once the growth on energy is mild at high energies in the present case. The 
agreement with the proton data is consistent and the extension to Ca and Pb targets is suitable 
since it is constrained by the DESY-HERA data and validity of the model in the energy range 
considered. Furthermore, the present investigation is complementary to those ones on heavy 
quark production [24] and nuclear structure functions [23] using the saturation model for nuclear 
targets. 

For light mesons production in 7p processes, the dipole Pomeron model predicts a larger 
growth with the energy than the saturation model due to the dominance of large qq pair sepa- 
rations in the saturation approach. This implies that Od% v oc cr , i.e. almost energy independent. 
It can be checked that the integration over dipole sizes in this dipole configuration gives an 
almost constant value, without logarithmic corrections as in DIS case. Differently, the dipole 
Pomeron model predicts a logarithmic dependence in the energy. In the nuclear case, this 
behavior implies a larger modification of the cross section in the dipole Pomeron model in com- 
parison with the saturation model. In particular, we have that for the nuclear exclusive vector 
meson photoproduction with A = Pb we predict that the difference between the results of the 
models should be a factor about 1.5. In contrast, for J/^> photoproduction, we have that the 
saturation model predicts a larger cross section for high energies. This is associated to the color 
transparency regime, present due to the small pair separation between charm and anti-charm. 
In this case we have a power-like behavior in contrast with a logarithmic dependence present 
in the dipole Pomeron model. 

A final comment on the small-t approximation considered here is in order. As shown in 
Ref. [49], the saturation effects play an important role in the t dependence of the scattering 
amplitude, mostly at large t. Therefore, the approximation presented here should be justified. 
Our master equation is exactly the same as in Refs. [49, 4], which reads, 




(33) 
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Figure 2: The total cross section for <fi and J/ty photoproduction on proton as well as for light 
(Ca) and heavy (Pb) nuclei. The solid lines stand for the QCD color dipole approach and 
the dashed ones for the soft dipole Pomeron approach. Experimental high energy data from 
DESY-HERA collider on proton target are also shown. 



where the squared momentum transfer is denoted by A 2 = —t. Basically, our expressions 
Eqs. (19) and (32) stand for the differential cross section at t — and further we have used 
the standard approximations for the small t behavior of the scattering amplitude on nucleon 
(proton) and nuclei. Namely, for the proton case one has considered the usual exponential 
parameterization which includes the meson slope parameter By. For the scattering on nuclei, 
this is accounted for by the nuclear form factor Fi(t), which includes the correct size of the 
nuclear target. Considering the saturation model, qq differential cross section is given by, 

(34) 



da q g (x,r,b) 
d 2 b 



1 - exp ( -- a dip (x, r) T(b) 



The remaining issue is what the accuracy of such an approximation concerning the sat- 
uration region for DIS with a nucleon/nucleus target. This can be addressed by looking at 
the Fourier transform of the differential dipole cross section for large dipole sizes (saturation 
limit). It was shown in Ref. [49] [see Figs. (16) and (17) in that paper], the saturation effects 
predict diffractive dips at large t. However, the pQCD picture remains quite the same at small 
t. Therefore, we believe that the integration on t of the complete expression, Eq. (33), should 
be not strongly sensitive to the large t region, once it is sub-dominant in the whole integrand. 
That is, we expect the effect on the total cross section is hidden in the integration, whereas it 
is important at large t in the differential cross section. 
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Figure 3: The total cross section for J/ty photoproduction for proton target (left panel) and for 
lead nucleus (right panel). The solid lines stand for the saturation model (SAT-MOD) and the 
dashed ones for the GLLMN model [50]. 

5 Summary and Conclusions 

In this paper we have calculated the nuclear cross sections for exclusive vector meson photo- 
production within QCD color dipole picture and Regge approach. These models predict cross 
section which does not violate the Froissart-Martin bound at high energies. Since they describe 
reasonably the experimental data for nucleon (proton) target, we are confident in extending 
these models for the nuclear photoproduction case. 

For the first approach, we have considered the saturation model, which is analytically simple 
and gives a good description of inclusive and diffractive ep experimental data. This model 
should be valid until the full non-linear evolution effects become important, which implies the 
consideration of the Pomeron loops beyond the multiple scattering on single nucleons estimated 
in the present framework. We have verified that the energy behavior is mild, mostly for J/ty 
where we would expect a hard behavior. This means that an important contribution also 
comes from large dipole configurations, related to soft domain. We predict absolute values for 
the cross section rather large, being about 4 mb and 0.12 mb for p and J/^f, respectively, for 
lead at W 7 jv ~ 1 TeV. These values are similar to those resulting from Ref. [9], except for 
the mild energy behavior for J/ty case presented here. Concerning the A-dependence, we have 
found a behavior proportional to A 2 / 3 (A 4 / 3 ) for the cross sections of light (heavy) mesons, in 
agreement with theoretical expectations associated with a transition to the black disc regime. 
In fact, for light mesons the nuclear shadowing in the scattering amplitude is stronger than a 
simple A 1 / 3 ansatz, as discussed in Ref. [23] when computing the nuclear structure functions. 
Namely, the saturation scale, which drives the A-dependence, between the proton and central 
nucleus is not simply oc A 1 / 3 , but has a prefactor which makes the result smaller. 

The results presented here can be contrasted, at least for the J/^f case, with the results 
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of Ref. [50] (hereafter GLLMN model), where the vector meson production (including DIS 
production) has been addressed in the color dipole picture. That analysis considers the Glauber 
approach and a numerical solution of the Balitsky-Kovchegov nonlinear evolution equation for 
the imaginary part of the dipole- nucleon scattering amplitude [51]. The comparison is shown 
in Fig. (3), where the solid lines represent the results from the saturation model (SAT-MOD) 
and the dashed lines are the numerical result of the GLLMN model. For the proton case (left 
panel), the saturation model gives a steeper growth on energy, whereas GLLMN produces a mild 
behavior at large energies. The upper/lower GLLMN curves stand for maximum and minimum 
values for the total cross section, obtained considering two different values for the correction 
factor Kp (we quote Ref. [50] for more details). It should be noticed that the behavior near 
threshold is also different in the two models. For the nucleus case (right panel), the behavior 
on energy remains basically the same as for the proton. The GLLMN model gives somewhat 
a slightly lower cross section at large energies, which it is about 8% smaller than the result 
for the saturation model at W 1 a ~ 400 GeV and presents a flatter behavior on energy. These 
features are directly associated with the different dipole-nucleon cross sections used in the two 
approaches. However, considering the relative errors of order of 15-20% in the predictions from 
Ref. [50], as stated by the authors, we can conclude that the agreement between the results is 
satisfactory, despite the distinct approximations made in the calculations. 

Concerning the nonperturbative approach, we have considered an unitarized Pomeron model 
and computed consistently its extension for nuclear targets. The energy dependence is loga- 
rithmic in any case and the nuclear effects seem to be stronger than in the saturation model 
for heavy nuclei. The nuclear dependence follows similar behavior as for the saturation model. 
However, it should be noticed we have used a quantum mechanical Glauber-Gribov calculation 
for nuclei targets in contrast with the results presented in Ref. [9]. As pointed out in Ref. [52], 
the present procedure gives a cross section higher than the classical mechanical model used for 
the predictions in Ref. [9]. For instance, it has been found in [52] a difference by a factor 2.5 
in p photoproduction at RHIC energies. 

Our results demonstrate that the experimental analyzes of nuclear exclusive vector meson 
photoproduction in the future electron-nucleus colliders eRHIC and HERA-A could be useful to 
discriminate between the different theoretical scenarios, mainly if heavy nuclei are considered. 
An alternative until these colliders become reality is the possibility of using ultraperipheral 
heavy ion collisions as a photonuclear collider and study vector meson production in this process. 
Moreover, such processes can be also studied outside the heavy ion mode. For instance, it has 
been discussed in Ref. [53] the photoproduction of heavy vector mesons in pp collisions at the 
Fermilab Tevatron and in the pp collisions at CERN LHC, since energetic protons also have 
large electromagnetic fields. These photoproduction reactions probe the gluon distribution in 
the proton at very small-x values [54], which open a new window to study parton saturation 
effects in exclusive processes. In a separated publication we will study these possibilities, 
considering the approaches discussed in this work. 
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